Monte Carlo study of exact S-matrix duality 
in non simply laced affine Toda theories 



M. Beccaria 
Dipartimento di Fisica dell'Universita and INFN 
Piazza Torricelli 2, 1-56100 Pisa, Italy 
e-mail: beccaria@hpth4 ■ difi. unipi. it 
(February 1, 2008) 

Abstract 

The (g^\d^) pair of non simply laced affine Toda theories is studied from 
the point of view of non perturbative duality. The classical spectrum of each 
member is composed of two massive scalar particles. The exact S-matrix 
prediction for the dual behaviour of the coupling dependent mass ratio is 
found to be in strong agreement with Monte Carlo data. 
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I. INTRODUCTION 



Affine Toda field theories are two-dimensional models described by the Euclidean action 
S = Jd 2 x (Id^ -d^+^r ±n a e^ (a) ^ . (1.1) 

The r-dimensional vectors {a' 4 '} are the simple roots of an affine Kac- Moody algebra 
and {n a } are positive integers depending on the algebra and satisfying 

^n a a a = 0, n = 1. (1.2) 

a 

The field is a set of r real scalar components. Finally, m and f3 are a mass scale parameter 
and the coupling constant. The Coxeter number is the positive integer h = n + ■ ■ ■ + n r . 

Under the transformation T : a — > 2a/\a\ 2 , the lattice of the simple roots transforms 
into the lattice of another affine algebra. The invariant algebras are called self-dual; they 
belong to the untwisted a-d-e series a£\ d£\ and to the twisted series a^. The other 
algebras are the pairs (bffl, a^n-i); ( c n \ ^2n-i)> (fl'lj df^), and (/j , eg 2 ' 1 ); they are invariant 
under T. 

At the classical level, affine Toda theories have no coupling; (3 can be scaled away, the 
spectrum is proportional to m, independent of j3 and moreover it is given by simple universal 
formulae in terms of the Coxeter number. The interest of the classical theory is that the field 
equations of motion admit a Lax pair and therefore there is an infinite hierachy of conserved 
currents with increasing spin. 

At the quantum level, this property is inherited in the form of a factorized S-matrix. 
The dependence on (3 which plays the role of Planck's constant becomes non trivial; on the 
other hand the parameter m becomes unphysical due to renormalization effects and only 
mass ratios are observables. 

Since the S-matrix is expected to be factorizable, its explicit form may be sought. One 
can make a guess and impose physical constraints like unitarity or crossing symmetry and 
the additional bootstrap principle. In the case of the self-dual theories, perturbation theory 
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suggests that the mass ratios do not renormalize. Indeed, the bootstrap equations close on 
an ansatz for the S-matrix based on the tree level spectrum and on the fusings allowed by 
the three-point couplings 

Perturbation theory and the structure of the bootstrap suggest conjectured expressions 
for the exact (3 dependence of the S-matrix which show a remarkable duality between weak 
and strong coupling in terms of the transformation f3 — > Air/ft. 

As pointed out in |4[], for the non self dual pairs the picture is more complicated. Mass 
ratios deform already at the lowest order of perturbation theory and the simplest ansatz for 
the S-matrix fails. 

However, a non trivial solution to the bootstrap equations can be found with the feature 
of predicting (3 dependent mass ratios [|5],|6|]. The predictions are then formally the same as 
in the classical theory, but in terms of a "renormalized" Coxeter number H((3). 

Again, the explicit non perturbative form of H(f3) is not known. The simplest conjec- 
ture |ID| , consistent with low order perturbation theory || and current algebra ]7| predicts 



a new kind of duality. Under (3 — > 47r/ ' (3 the S-matrices of the pair members get exchanged. 
Hence, the strong coupling regime in one theory should be given by the weak coupling regime 
in the other. 

In [H a Monte Carlo study of duality in the pair (g^\ d^) was performed by mean of the 
Metropolis algorithm. The authors determined the mass ratio in the theory over a wide 
range of couplings and they did find agreement with the duality conjecture. Specifically, 
they checked that the mass ratio in ranged between its classical values and the classical 
value of . 

In this paper I carried over the above simulation on larger lattices with higher statistics 
in order to pin down the precise dependence on (3. Moreover, I have used the Hybrid Monte 
Carlo algorithm ||13|. Finally, I have extended the simulation to the d\ theory in order to 



have a complete picture. 

The plan of the paper is the following: Section [TI| describes the pair 

(^ 1} ,4 3) ) ; Section Q 



shows the one loop deformations of the mass ratios; Section IV gives some detail on the 



numerical simulation; finally, in Section [V] the results are discussed. 



II. THE DUAL PAIR (G { 2 1} , D ( f) 



The pair (<7 2 , d± ) has r = 2 and its action is 

f fl m 2 2 1 

S = J d 2 x^ -<9 M ■ <9 M + — ^ n a exp(/3 a a • 0) j , 



(2.1) 



The integers n a and the simple roots are 
gP :n = {2,3,1}, a 

df : n = {2, 1, 1}, a 



{(V2,0),(- 



(2.2) 
(2.3) 



The two sets of roots are related by the duality a — > 2a/\a\ 2 . 

The two corresponding models are very different at the tree level. The explicit expansion 
of the mass-potential exponential term up to the fourth order for the gtp model is 



V(0i,0 2 )=m 2 (3 2 + 2 ) + 

9 03_ 90i0 2_ 2v /3 3 

9^ 

| m2/?2 270f + 18 2 2 + 8 y^fc 0| + 70| + 0(/33) 

72 



(2.4) 



(3) 

For the c?4 model we utilize the tree level mass eigenstates by transforming the fields 



(f)^R<f>, R 
and obtain the expansion 



^ cos 6 sin 6 ^ 



sin# cos 9 



5tt 
12 



(2.5) 



V(<f>!, 2 ) = m 2 ((3 - V3) 2 + (3 + \/3) 2 ) + 

+m 2^ _ + Mfo _ y^M^g + + v^^i^i + 0| + VM) + ( 2 6 ) 



2 ,2 / ' rl 5 <j)\ 0? fa 0? 2 , 3 2 2 01 0| 01 0| 



\ 8 4^ 2 ^3 4 



, 70| | 50^ 



8 4-^3/ 



+ 0(/3 3 ) 
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As one can see, the sets of possible fusings are completely different and duality is far from 
being obvious. The classical mass ratios are 



m 2 
nil 



m 2 
mi 



,(3) 

H 



V3_ 
\V3 



(2.7) 



1 V2 

which agree with the general formula in terms of the Coxeter number h (6 for , 12 for 



4 3) ) 



m 2 sin(27r//i) 
nil sin(7r//i) 



= 2 cos(n/h) 



(2.8) 



The duality conjecture states that the correct quantum ratio gtp is given by substituting 
h -> H(P) in the model g { 2 l) and h -> H(Att/(3) in the df model. The form of H((3) is 
constrained but not fixed by perturbation theory and the conjectured expression is 



ff( 0= 8+ l +w 

Let us clarify these statements by considering the one loop mass ratios. 



(2.9) 



III. ONE LOOP MASS RATIOS 

Let us denote the three diagrams of Figures (I-II-III) by 

r (l) p(2) p(3) 

1 abci 1 abcdi 1 abed \°- L ) 

where a, b, c and d are particle labels in the range {1,2}. The mass ratio is observable 
since renormalization amounts to a normal ordering of the exponentials and its effect is a 
redefinition of the bare mass. We must check that in a bare renormalization scheme all 
the divergent tadpole graphs cancel. Let us utilize dimensional regularization and let us 
introduce 

Zl = l (270V +m 2- (3 ' 2) 

The pole part of Zi is mass independent, hence the cancellation is a matter of couplings. At 
the one loop level the mixed propagator corrections are irrelevant and we can restrict to the 
diagonal ones. Let us write the interaction lagrangian in the form 
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V(<f> u <h) = \{m\(t>l + ml4>l) + VukPI + V ll2 4>\<p 2 + ■■■. (3.3) 
Then, the corrections to the propagator of particle 1 are 



1 in 


= -12 Vim £i, 




1 112 


= — 2 V1122 Z 2 , 




F (2) 
1 1111 


= 18 V? u Z x ml 


2 


F (2) 
1 1112 


— 6 Vin V122 ^2 




F (2) 
1 1121 


= 2 lf 12 Zi m 2 " 2 




F (2) 
1 1122 


= 6 V112 V222 -^2 


m 2 " 2 



(3.4) 



The corrections to the propagator of particle 2 are obtained by exchanging the 1 and 2 
labels. If we denote the full divergent correction by 

X m 2 _ p(l) , p(l) , p(2) p(2) p(2) p(2) ,„ ^ 

— 1 111 + 1 112 t 1 mi -|- 1 1112 r 1 1121 1 r 1122 

then the desired cancellation is equivalent to the condition 

5m\ bm\ 



1 1 



(3.6) 



which is indeed satisfied by the couplings of the two theories which can be read in the 
expansions of the previous section. 

Besides the consistency check, let us turn to the mass ratio deformation. At one loop, 
we must determine the quantity 



^ = ^R-^ • (3.7) 
m 2 m 2 y mf m 2 ) 

Let us introduce the finite integral 

2 r d 2 q 1 
lAP ] = J (27r) 2 (g 2 + m 2 )((g + p) 2 + m 2 )- { } 

Then the finite contributions to the propagators of particle 1 are 
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rfAi = isv&ZuCp 2 ) 

^1122 = 21^22^22 (p 2 ). 



(3.9) 



Evaluation on the tree mass shell gives 



bm\ = 18V? u Z u (-ml) + 4lf 12 Z 12 (-m 2 ) + 2lf 22 Z 22 (-m 2 ) 



(3.10) 



with analogous expressions for the particle 2. We need only the following particular values 

1 1 



Zn(—m, 



2 ; 



27r A /m 2 (m 2 — 4m 2 ) 



=ArcTanh 



\ 



m 2 — 4m? 
m 2 



(3.11) 
(3.12) 



and the final result is 



^ = _^ + 0(/? 4 ); d. 



1 



(3) . 



m| 12^3 
The renormalized Coxeter number is thus 



<7 2 ': #(/?) = 6 + ^ + 0(/3 4 ); 4" 



and a consistent, simple and natural conjecture is 



H{(5) = 12 - -f- + 0(/3 4 ) (3.14) 



#2 



#(/?) = i2o(/3); 4 



= Fo(4tt//?) 



(3.15) 



where 



H Q {/3) = 6 + 



6<H n < 12. 



(3.16) 



1 + /? 2 /12tt 

The result for is in agreement with that quoted by |J . The result for g^ 3 " 1 gives perturba- 
tive support to the duality conjecture f3 — > in/ft. We remark that the discrepancy with || 
is due to the fact that they use the form of H Q which is correct for the simply laced models. 
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IV. DETAILS OF THE SIMULATION 

The lattice action for the pair (g^\ df^) expressed in terms of pure numbers is 

f 1 m 2 3 1 

^Toda= \t; Y ( ( Pn+^~ <Pn) 2 + ~wJ2n a exp(P a a - , 0=(0 1 ,0 2 ). (4.1) 

nesites I Z /u=l,2 P a=l I 



I have simulated the Toda theory with the Hybrid Monte Carlo algorithm (see [13 for the 
details). Let us consider the extended action 

S = S p + Sroda, (4.2) 

S P = \ ■ VT n , 7T = (7Ti,7r 2 ). (4.3) 

The free parameter of the algorithm are Nh mc and e. The first is the number of molecular 
dynamics steps. The second is the time step in the integration of the equations of motion 



(4.4) 



2 

m 



= XX^W ~ 2< ^« + ~ -rJ2n a a a exp(pa a ■ <p n ). (4.5) 

fj, Pa 

The vacuum expectation value of the field is a non physical quantity. However, it is inter- 
esting to measure it since it is an indicator of thermalization and also because it is in a sense 
a dynamic minimum of the Toda potential. 

Mass ratios can be determined by studying the eigenvalues of the two point function 

(0\M0)®Ar)\0) - mmiOl^lO) (4.6) 

where t is the lattice time ranging from to T and the wall field $«(£) is obtained by 
averaging over space. 



V. RESULTS 

I have used a 80 2 lattice for all /3s because the correlation lenght may be adjusted by 
varying m. The continuum mass ratio is independent of the bare mass m. However, on a 
finite lattice, it must be chosen in order to have correlation lenghts large with respect to one 
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lattice spacing and small compared to the lattice size. This is the correct procedure which 
minimizes discretization and finite size corrections. Thanks to the work of || I had good 
values in the case of the theory. In the other model, I started with the same values of 
m adjusting them for some f3. 

I utilized different measurement of the wall-wall two-point function for each bin in the 
separation r. This is necessary in order to avoid strong correlation between data. 

Table | shows the Hybrid Monte Carlo parameters which we found to be optimal for 
each couple ((3, m). The time step e must be reduced almost exponentially as (3 is increased. 
This is reasonable since at larger (3 the potential profile becomes steeper. 

Table H shows the measure of (O|0j|O) which can be useful as a check of the code and 
which is needed in order to subtract the two-point function. 

Table |T| shows the two lattice masses, their ratio and the conjectured prediction. 

Finally, tables [TV], [V|, |\T] show the same results in the case of the d\ model. 

Figures I-II-III show the self energy diagrams which are needed in order to compute the 
one loop mass ratio deformations. 

Figure IV shows a summary plot of the measured mass ratios in the two models together 
with the conjectured ones and the asymptotic values holding in the classical limit. 



VI. CONCLUSIONS 

In this paper I have investigated numerically the conjectured duality in the pair (g^\ df ) 
of non simply laced affine Toda theories. I have shown that the (3 dependence of the mass 
ratios in does follows the behaviour conjectured in || and that the data of d^ agree 
with the (3 — > An / (3 duality. 

As in the case of more realistic field theories like QCD, the numerical approach could 
be useful in studying other interesting features of quantum Toda theories. For instance, 
one could try to find direct evidence of the boundary bound states which appears when the 



theory is restricted to a half-line flIT|]; work is in progress on this topic. Moreover, it could be 
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valuable a non perturbative study of the solitons which appear at imaginary /3, and which 
suggests that a unitary theory can ultimately be found by restricting the state space of the 
hamiltonian [113]; their stability is indeed still questionable ||14j|. 
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CAPTIONS 

Fig. I-II-III : self energy one loop diagrams. 
Fig. IV : summary of the numerical results. 
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TABLES 
TABLE I. g£>: HMC parameters. 



P N hmc e 


1.0 




10 




0.1 


2.0 




10 




0.08 


3.5 




10 




0.07 


5.0 




10 




0.06 


10.0 




5 




0.05 


20.0 




10 




0.025 




TABLE II. 


gip bare mass and dynamical minimum. 




P 


m 






<02> 


1.0 


0.1 




-0.0982(6) 


0.225(1) 


2.0 


0.1 




-0.1542(4) 


0.3768(8) 


3.5 


0.05 




-0.2073(5) 


0.5300(9) 


5.0 


0.01 




-0.2752(4) 


0.6866(8) 


10.0 


5E-5 




-0.3396(7) 


0.843(1) 


20.0 


5E-7 




-0.2362(4) 


0.7023(7) 
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TABLE III. g^: mass ratio. 






m\ 


m 2 


R 


R* 


1.0 


0.1595(2) 


0.27839(5) 


1.745(2) 


1.74509 


2.0 


0.2133(2) 


0.3784(1) 


1.775(2) 


1.77605 


3.5 


0.2262(5) 


0.41192(5) 


1.821(5) 


1.82579 


5.0 


0.1615(4) 


0.3004(1) 


1.861(5) 


1.86150 


10.0 


0.1379(3) 


0.2630(2) 


1.908(5) 


1.90870 


20.0 


0.2896(3) 


0.5576(1) 


1.926(3) 


1.92562 


TABLE IV. 4 3); HMC parameters. 







Nhmc 




e 


1.0 




10 




0.07 


2.0 




10 




0.08 


3.5 




10 




0.05 


3.5 




10 




0.06 


5.0 




10 




0.05 


10.0 




10 




0.02 


10.0 




10 




0.03 


20.0 




10 




0.0125 
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TABLE V. d£> bare mass and dynamical minimum. 






m 




{<h) 


1.0 


0.1 


0.1530(5) 


-0.1930(9) 


2.0 


0.1 


0.2014(3) 


-0.2210(6) 


3.5 


0.05 


0.2305(2) 


-0.2174(4) 


3.5 


0.01 


0.3163(5) 


-0.3687(9) 


5.0 


0.01 


0.2701(3) 


-0.2552(5) 


10.0 


5E-5 


0.3061(3) 


-0.2729(6) 


10.0 


1E-7 


0.4410(8) 


-0.504(1) 


20.0 


5E-7 


0.2488(2) 


-0.1489(4) 



TABLE VI. 4 : mass ratio. 






mi 


m 2 


R 


R* 


1.0 


0.2166(1) 


0.4149(1) 


1.916(2) 


1.91665 


2.0 


0.3363(2) 


0.6323(1) 


1.880(1) 


1.88120 


3.5 


0.3427(1) 


0.6280(1) 


1.832(1) 


1.82840 


3.5 


0.2150(1) 


0.3931(3) 


1.829(1) 


1.82840 


5.0 


0.3850(2) 


0.6903(4) 


1.793(2) 


1.79398 


10.0 


0.3506(6) 


0.6155(3) 


1.756(4) 


1.75193 


10.0 


0.1201(3) 


0.2106(1) 


1.753(6) 


1.75193 


20.0 


0.5727(6) 


0.9966(6) 


1.740(2) 


1.73740 
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a 




1.95 



1.90 



1.85 



1.80 



1.75 




1.70 1 L — ^ ^ 

0.0 5.0 10.0 15.0 20.0 

P 



